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We study the exclusive semileptonic decays B s — > D* £i) and B s — > D* x £v, where 
p-wave excited D* and states are identified with the newly observed D s j(23l7) 
and -D s j(2460) states. Within the framework of HQET the Isgur-Wise functions 
up to the subleading order of the heavy quark expansion are calculated by QCD 
sum rules. The decay rates and branching ratios are computed with the inclusion 
of the order of 1/mg corrections. We point out that the investigation of the B s 
semileptonic decays to excited D s mesons may provide some information about the 
nature of the new D* sJ mesons. 
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I. INTRODUCTION 

Recently BaBar collaboration reported a narrow state with J p = + at a rather low 
mass 2317 MeV decaying to D s 7r° This observation was quickly confirmed by CLEO, 
and another narrow state with J p = 1+ at 2457 MeV was found in the D* s n° channel 0. 
Both states have been confirmed by Belle |3J]. Since these states lie below DK or D*K 
threshold, they decay instead via isospin-violating transitions. 

These newly observed states have attracted much attention because their measured 
masses and widths do not match the predictions from potential-based quark models 0|. 
To resolve the discrepancy, many theoretical speculations have appeared in the literature. 
Bardeen et al. interpreted them as positive-parity cs (0 + , 1 + ) spin doublet of the Eh and D* 
negative-parity cs (0~, 1~) ground states in the framework of chiral symmetry [5J]. Based 
on the quark-antiquark picture, various theoretical models are modified to accommodate 
the low masses and the narrow widths for the new states QCD sum rule 

analysis in supports the quark-antiquark postulation, both new states, -D s j(2317) and 
-D s j(2460), are identified to be the cs excited + and 1 + states in the ji = \ + doublet. 
Apart from the quark-antiquark interpretation, the _D*j(2317) meson has been interpreted 
as a DK molecule a D s i{ molecule 0^ four-quark state [lij], and a mixing of the 
conventional state and the four-quark state |l5l |. 

Motivated by the interpretation that the two new states are the excited p-wave D s 
states belonging to the (0 + , 1 + ) doublet with ji = i + in jll|, it is worthwhile to investigate 
the B s exclusive semileptonic decays to this doublet, ie. B s — > D* £v and B s — > D* sl Ei>. 
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The search of these decay modes may provide some information to understand the nature 
of these new states and to clarify the controversy. Although current B factories do not 
produce B s meson, these decays can be studied at future hadron B factories. 



In this paper we shall use QCD sum rules 16, ^ in the framework of the heavy quark 
effective theory (HQET) |19| to study exclusive semileptonic B s decays to excited spin 
symmetry doublet (D* Q , D* x ) mesons. HQET is a useful tool to describe the spectroscopy 
and weak decays of hadrons containing a single heavy quark, it provides a systematic 
method to compute the properties of heavy hadrons via the 1/rriQ expansion, where niq is 
the heavy quark mass. The study for the non-strange B semileptonic decays into charmed 
meson doublet (0 + ,l + ) in HQET can be found in the literatures by using various ap- 
proaches, including HQET-based considerations [20 , EH , QCD sum rules j^J EE E^J and 
various quark models [25, 2f], 27, 28, 2^E3|- m this work, we shall calculate the weak decay 
elements for B s — > D* , D* sl in HQET. Considering that most of the phase space for these 
decays is near zero recoil, we shall include the Kqcd/^q corrections in the application. 

The paper is organized as follows. In Sec. |H]we review the formulas for the matrix ele- 
ments of the weak currents including the structure of the ^QCD/ m Q corrections in HQET. 
The QCD sum rule analysis for the leading and sub leading Isgur-Wise functions is pre- 
sented in Sec. IIHI Section HY1 is devoted to numerical analysis and the applications to 
decay widths. This section also includes a brief summary. 



II. DECAY MATRIX ELEMENTS AND THE HEAVY QUARK EXPANSION 



The matrix elements of vector and axial vector currents (V^ = cj^b and = 07^75 b) 
between B s meson and excited D* or D* sl mesons can be parameterized as 

(D* s0 (v')\V»\B s (v)) = 0, 

(D* s0 (v')\A»\B s (v)) = g+W + v'^+g-W-v*), 
(D* sl (v',e)\V»\B s (v)) = g Vl e*» + (g V2 v» + g v ^) e>V , 



(D* sl (v',e)\A»\B s (v)) = ig A er*'&ri 



(1) 



Here form factors g^ are functions of the dot-product, y = v ■ v', of the initial and final 
meson four-velocities. The differential decay rates expressed in terms of the form factors 
are given by 



dr 



D* 
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± rl (y 2 - If' 2 (1 + r )g + - (1 - r Q )g. 
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(2a) 



B 



dy 



487T 3 

+2(l-2r iy + rl)[g Vl + (y 2 -l)g 2 A 
where r = m D * /m Ba and n = m D * lm B . 

sO ' s si ' 3 



(y - g Vl + (y 2 - i)(sv 3 + r i9v 2 ) 



(2b) 
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The form factors gi can be expressed by a set of Isgur-Wise functions at each order 
in Aqcd/^c> This is achieved by evaluating the matrix elements of the effective current 
operators arising from the HQET expansion of the weak currents. This problem has been 
discussed previously in Ref. [21[, we outline here the analysis for the sake of completeness. 

One introduces the matrix representations 



H v 
K v 



i + i> 

2 



(3) 



composed from the fields P v , P*^ and P' vl P'j 1 that destroy mesons in the doublets j[ = | 
and | + with four-velocity v in HQET, respectively. At leading order of the heavy quark 
expansion the hadronic matrix elements of weak current between the states annihilated by 
the fields in H v and K v > are written as 



h^Th^ = ((w)Tr{K v ,TH v y 



(4) 



where is the heavy quark field in the effective theory and ( is a universal Isgur-Wise 
function of y. 

At order Aqcd/^q there are contributions to the decay matrix elements originating 
from corrections to the HQET Lagrangian 



5C 



2m Q L 



nW) + nW) 



(5) 



and originating from the matching of the b —> c flavor changing current onto those in the 
effective theory 



cTb = T$ I r 



(6) 



where D is the covariant derivative. The matrix elements of the later operators can be 
parameterized as 



h$iD x Th® = Tr \s{ c) K v ,TH v 



h^riDxh® = Tr{<Sf K V ,TH V }. 



(7) 



The most general decomposition for is 



S 



(0) 



ti Q) v x + cPv x + cf ^ . 



The functions Q depend on y and have mass dimension one. The translation invariance, 
id v (h^/Thv ) = (Atv — A.'v' u )h^/ Thv , and the motion equation for the heavy quark, 
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iv ■ Dh^ = 0, result in the relations between the form factors (- w> |21l 



(Q) 



C ( C ) + C? ) 



-(c) 



AC, 
yA'-A 

y + i 



d c) +d b) 



C-Ci 



Ac) 

S3 



-A'C, 
_ 2/A' - A 
y+1 



d c) + d 6) 



C-(y-i)Ci 



o 



(9) 



where A (A') = mui+nM') — ttiq is the difference between heavy ground state j\ 



(excited j[ 



meson and heavy quark masses in the rriQ — > oo limit. These relations 



show that all corrections to the form factors coming from the matching of the weak currents 

(c) 

in QCD onto those in the effective theory are expressible in terms of Q and Q . 

The matrix elements of Aqcd/^q corrections from the insertions of the kinetic energy 
operator Okin have the structure 



d 4 xT O 



kin,i)' 



i J d 4 xT j 



(o) 
(o) 



xtTi \R v ,TH v \. 



(10) 



The functions Xkeiv) have mass dimension and effectively correct the leading order Isgur- 
Wise function ((y) since the kinetic energy operator does not violate heavy quark spin 
symmetry. 

There are Aqcd/wq corrections associated with the insertion of chromomagnetic op- 
erator O mag . The QCD sum rule approach for the semileptonic B decays to ground state 
and excited D mesons shows that the functions parameterizing the time-ordered products 
of the chromomagnetic term in the HQET Lagrangian with the leading order currents are 
neglig ibly small 0,0- This is in agreement with the results obtained from the HQET- 
motivated considerations [ijj and relativistic quark model [EJ . Therefore, we shall neglect 
the chromomagnetic correction hereafter. 

Summing up all the contributions the resulting structure of the decay form factors is 



9+ 
9- 
9a 

9Vi 

9v 2 
9v 3 



2(y-l)Ci-3C 



3/A' - A 
y + 1 , 



A'(2y + 1) - A(y + 2) 



C + £c Xke + Zb Xke 



yA'-A 



C + xL 



M(2y + 1) 



y+1 

■My + 2) 



C-2(y-l)Ci 



y + 1 



C-2(y-l)Ci-xL 



(y - 1) C + e c [(yA' - A)C + (y- i)xL 

-e b { [A'(2y + 1) - A(y + 2)}(- 2(y 2 - 1) d - (y - l)xL} 



2e c ( u 
-C-e 



yA'-A 



C + 2Cl + Xke 



y+1 

A'(2y + l)-A(y + 2 
y + 1 



+ 



C-2(y-l)Ci-Xko 
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where Eq = l/(2mg) and the superscript on Q c is dropped. In the following sections we 
shall employ the QCD sum rule approach to calculate the leading and subleading Isgur- 
Wise functions. 



III. FORM FACTORS FROM QCD SUM RULES 
A. Leading Isgur-Wise function £ 



The proper interpolating current J^p^ 3 for a heavy mesonic state with the quantum 



number j, P, j% in HQET was given in 32j. These currents were proved to satisfy the 
following conditions 

m^mf.P'Ji) = fPn^PP'^a'f 1 ^ , (12) 

/(()7-(./;;;r;'xr)./;/;,;J(())) |o> = 5 jf 8 PP ,5 jefl (-iys g ^-.-g^ 

x J dt5(x - vt) U Pth (x) (13) 

in the uiq — > oo limit. Where 7] ai " a i is the polarization tensor for the spin j state, 
g" 13 = g alS — v a v lS is the transverse metric tensor, S denotes symmetrizing the indices and 
subtracting the trace terms separately in the sets (a\ ■ ■ - ay) and ■ • ■ Pj), fpj e and Hpj e 
are a constant and a function of x respectively which depend only on P and jg. 
The local interpolating current for creating 0~ pseudoscalar B s meson is taken as 

4,-1/2 = \J\ hvlss , ( 14 ) 
and the local interpolating currents for creating + and 1 + (D* Q , D*^ mesons are taken as 



1 

3 



J o,+,i/2 = ^ h v{-i)PtS 



4,+,i/ 2 = 4= K^ti-i) p t s , (15) 

In order to calculate this form factor using QCD sum rules, we follows the same proce- 
dure as [33| to study the analytic properties of the three-point correlators 



i 2 J d 4 xd 4 ze^ k '- x - kz \0\T(j ^ 1/2 (x)Jji(0) Jl_ >1/2 (z))\0) = 5(w, u/, y)C^ , (16a) 
i 2 / d 4 xd 4 ze^- k - z \0\T(j^ y2 (x) J^JO) 4 t _ tl/2 (z))\0) = S(«, u/, y)C^ A , (16b) 



where Jy = h^v')^^ h(v) and J A = h{v')^^ r )^, h(v) are leading order vector and axial vector 
currents, respectively. The variables k, k! denote residual "off-shell" momenta which are 
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related to the momenta p of the heavy quark in the initial state and p' in the final state by 
k = p — ttiqv, k' = p' — rriQ'v', respectively The Lorentz structures, Cy,A, have the forms 

L\ = v» - v" , C% = -ie^v a v'p , C% = (y- l)gjT ~ t/X , 

where vf = v@ = v a — yv' a . 

The coefficient S(ou, uj', y) in (fTBj) is an analytic scalar function in the "off-shell energies" 
uj — 2v ■ k and uj' = 2v' ■ k! with discontinuities for positive values of these variables. 
By saturating the double dispersion integrals for the correlators in (jlfij) with physical 
intermediate states in HQET, one finds the hadronic representation of the correlator as 
following 

n hadro (u, u,y) = —j . Wo t. ; r-r + higher resonances . (17) 

(2A — uj — le) (2 A' — uj ' — it) 

As the result of equation (fT2j) . only one state with j p = 1 + contributes to ()17j). the other 
resonance with the same quantum number j p and different ji does not contribute. 

On the other hand, the correlator can be calculated in QCD in the Euclidean region, 
i.e., for large negative values of uj and uj', in terms of perturbative and nonperturbative 
contributions. Furthermore, the nonperturbative effects are able to be encoded in vacuum 
expectation values of local operators, the condensates. Hence one has 

U(uj,uj',y) = / dvdv — - — + n con d + subtractions . (18) 

J [v — uj — ie){u' — uj' — ie) 

The QCD sum rule is obtained by equating the phenomenological and theoretical ex- 
pressions for S. In doing this the quark-hadron duality needs to be assumed to model the 
contributions of higher resonance part of Eq. (JTTj). Generally speaking, the duality is to 
simulate the resonance contribution by the perturbative part above some threshold ener- 
gies. In the QCD sum rule analysis for B semileptonic decays into ground state D mesons, 
it is argued by Blok and Shifman in ji^] that the perturbative and the hadronic spectral 
densities can not be locally dual to each other, the necessary way to restore duality is to 
integrate the spectral densities over the "off-diagonal" variable uj_ = [y — v')/2, keeping 
the "diagonal" variable uj + — {y + v')/2 fixed. It is in uj + that the quark-hadron duality 
is assumed for the integrated spectral densities. We shall use the same prescription in our 
application. 

In order to suppress the contributions of higher resonance states a double Borel trans- 
formation in uj and uj' is performed to both sides of the sum rule, which introduces two 
Borel parameters T\ and T2. For simplicity we shall take the two Borel parameters equal: 
T 1 = T 2 = 2T. 

The calculations of p peTt and II con( j in HQET are straightforward. In doing this, for 
simplicity, the residual momentum k is chosen to be parallel to v such that k^ = (k ■ v 
(and similar for k') in the theoretical calculation. Since we deal with Qs states, the light 
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quark mass shall be included in calculating the spectral function and condensates. For the 
condensates we confine us to the operators with dimension D < 5 in OPE. By performing 

;he Taylor expansion around x p = 0, we obtain the expansion for the quark-condensate 

3? 



^j a {x)D a tfj p (y) :) = -|^-[im g (7 a ) pCT - 5 pa m 2 q (x Q - y a )\ + q9s(7 ^ Gq {(g,^ _ ^a afl ) pa 



(x" - y^) + — m q [^ a ) pu {x - yf + 2($6 - ^) P<J (x Q - y a )] 

+2(x a - y a )y p x v )] } . (19) 

After making double Borel transformations in the variables u and u' and changing the 

integral variables co> + = (z/ + i / ')/2, u;_ = ( ) — z^ / )/2, one obtains the sum rule for 

£ as follows 



C(y) /-i / + ,| e" (A+A ' )/T = ^7-^ / ^ + e~- /T [w» + 3(1 + y)(m. w J + 



m 2 ^ 



(ss), 2 2/ + !\ 1 2/-\ 1 + 2//i 5m s m 2 15y + 7, 

Jm, — m Q mn(ss) 1 — 

g V s s rp ) 12 o\ / T \ 8T 48 T 2 ' 

-GG) —— + —-(—GG) [2 lE - In — + In — — ) , 20 

7T V + 1 32T 7T u 2 2 7 



192 x vr ' y + 1 32T it ' v /i 
where m 2 , (ss) = (sga pu G^ w s) with m 2 , = 0.8GeV 2 . 

B. Subleading Isgur-Wise function 

In order to derive the QCD sum rule for the subleading form factor (i(y) defined in 
we consider the following three-point correlation functions 



Sg^wy.y) = i 2 J d 4 xd*ze« M ^ (0|T (j , + , 1/2 (x) h$ ip^ l5 /if (0) j£_ fl/a (z)) |0> 
= H (o;,a; ,y)£^ , (21a) 

H^(u,,u/,</) = * 2 / Jxd'zeW*-'-') (0\T (r w/2 (x) h { :?ipr l5 h^(0)jl_ >1/2 (z)) |0) 
= S^u/^KT, (21b) 

s&(a;V,v) = * 2 / «^ ,(t '' 1 '- k) (o|T (jr, +il/2 (x) htfip-fWio) Jj f _, x/a (*)) |o) 

= S 2 (w, w', y)£^y + S 3 (w, u', y)£% v , (21c) 
where the Lorentz structures have the forms 

au = «" + . ^ = sTG/ + 1) - «V , £gv = «x - « V . 
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The coefficient functions Ei(u>,uj',y) can be expressed in terms of perturbative and non- 
perturbative contributions in QCD theoretical calculation. These functions are used to 
construct the sum rules needed. 

By saturating the double dispersion integral for the three-point functions in (}2~T]) with 
hadron states and using Eqs. ©-© and (fT2"j) . one can isolate the contributions from the 
double pole at uj = 2 A, uj' = 2A': 

/ 1 f 12 

= (2A-o;-^)(2A' 2 -a;'-i e )[ 2(y - 1)Cl(y) - 

-3^^C(y)] + ■ ■ ■ , (22a) 

EZ(w,v',y) = f ~\ U T r yA '~ A C(y)C7 + --- , (22b) 

1AK ' ,yj (2A - uj - ie)(2A' - uj' - ie) y + 1 SVy; A V ' 

i \ /-,|/+,i/2 yA' — A 

(2A — uj — xe) \2A' — uj — te) y + 1 s 

-2d(y)£gV] + --- ■ (22c) 

From (|2*T|) and (|2*2*|) one can see that in the case of S 3 , the residue of the pole is proportional 
to the universal function Ci(y)> the pole contribution to Hi and S 2 is related to ((y)- While 
for S the residue of the pole is proportional to both ((y) and Ci(y)- QCD sum rule is 
obtained by equating the phenomenological and theoretical expressions for 5. Therefore, 
the sum rule for the subleading form factor (i(y) can be constructed either from S3 or from 
So- One can also yield the sum rule for the leading form factor in the form (yA' — A)((y) 
from S , Hi and S 2 , respectively. 

We shall focus on the coefficient function E 3 (uj,uj',y) to construct the sum rule for the 
subleading form factors (i(y). One obtains the sum rule for (i(y) as follows 



Ci(v) /-,*/+,* e - (A+A ' )/T = "1^(^1)3 y 2ms ^ + e-^K + 4m s ( 2/ + l)(2 2/+ l) 2 4] 

m^M _ T . (23) 

96 T 96 n 1 (y + l) 2 y ' 

Moreover, the sum rule for the combination (yA' — A)((y) can be obtained independently 
from the coefficient function Si and S 2 in (|2"T|) together with (|22|1. respectively. We have 
double checked that the resulted sum rule for ((y) has the same form as (|20|1 . The above 
consistency checks confirm that our method is consistent with the general analysis of Ref. 
jijj described in Sec. HU 



C. QCD sum rules for 77? 



For the determination of the form factor 77^, which relates to the insertion of Aqcd/^c 
kinetic operator of the HQET Lagrangian, one studies the analytic properties of the three- 
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point correlators 

i 2 J d*!^^ (0|T (j 0l+ ,i/ 2 (x') O^(z)^(0) Jj,_, 1/a (x)) |0> 

= £(u;,u/,y)££, (24a) 
* 2 | ^sdV^ze^-**) (0\T (r w/2 (x') 0^ v ,(z) J» A (0) Jl_ A/2 (x)) |0) 

= S(w,o/,y)££ 1 . (24b) 

By saturating (|2*4^) with physical intermediate states in HQET, one can isolate the 
contribution of interest as the one having poles at to — 2A, to' = 2A'. Notice that the 
insertions of the kinetic operator not only renormalize the leading Isgur-Wise function, 
but also the meson coupling constants and the physical masses of the heavy mesons which 
define the position of the poles. The correct hadronic representation of the correlator is 

- / / \ /-,l/2/+,l/2 I c f \ , 

n hadro (o;,u;,y) = j-= — -j- — [r] ke {y + 

(2A — uj — te)(2A' — to' — xe) \ 

{G+,1/2 + r) t ~r~~~T~ ■ ) C(!/) J + higher resonance , (25) 
7 2A' — uj' — le J 



where Kpj e and Crp f are defined by 



32 



3 

W) 



V,PJi\OSZv\3,P,3i) = K p,u> 

/ ^ ^i:? ; :,('')-/;:;;r ! «'!. /./'../,} = /p A c?^-^ . (26) 

Within the same procedure one finds the sum rule for r] ke as 



bUv) + (Gf ll/2 + ^)C(y)] /-,i/ 2 / + ,i/2 e" (A+A ' )/T 



^4 

2T 

3 3y + 2 r ... .^/rr.i , 4. 
16vr 2 (y + l) s 



5 2 (ss) l.a s . 3y + l 
+ 32 + 96 { ~ GG) (^Tiy T ' (27) 

From the consideration of symmetry, the sum rule for 7] ke that originates from the 
insertion of AqcdA 7 ^ kinetic operator of the HQET Lagrangian is of the same form as 
in (|2*7|). but with the HQET parameters G+i/ 2 and if+,1/2 replaced by G^ l , 2 and K^ t y 2 , 
respectively. The definitions of G K : 1 , 2 and i^-,1/2 can be found in Eq. (J2EJ). 

We end this subsection by noting that the QCD 0(a s ) corrections have not been included 
in the sum rule calculations. However, the Isgur-Wise functions obtained from the QCD 
sum rule actually are a ratio of the three-point correlator to the two-point correlator 
results. While both of these correlators subject to large perturbative QCD corrections, 
the remaining corrections to the form factors themselves may in fact be small because of 
cancellation. This is what happens in the case of analysis for B semileptonic decay to 
ground state and excited charmed mesons [^E^]. 
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IV. NUMERICAL RESULTS AND PREDICTIONS FOR SEMILEPTONIC 

DECAY WIDTHS 

We now turn to the numerical evaluation of these sum rules and the phenomenological 
implications. In order to obtain information for £(?/), Ci{v) and Vkeiv) fr° m the sum rules, 
we need the related decay constants, /-,i/2 and /+,i/2, defined in (|T2j) as input. The QCD 



sum rule calculations for the correlators of two heavy-light currents give [111 l38j|: 



+^<-> (l-3T + 3#J - T6P ( V GG > (** " 1 " ' V 1 ' <28) 



64vr J2m 3 



(5 a >(l - ^ + ^) + ^ («> - ^ (« S G 2 ) . (29) 

These two-point sum rules can be used to eliminate the explicit dependence of three-point 
sum rules on /_ 1/2 and /+,i/2, as well as on A and A'. This procedure may help to reduce 
the uncertainties in the calculation. 

For the QCD parameters entering the theoretical expressions, we take the standard 
values 



(qq) = -(0.24±0.01) 3 GeV 3 , 

GG) = (0.012 ± 0.004) GeV 4 . (30) 



The strange/nonstrange condensate ratio is adopted as (ss) = (0.8 ±0.1) (99). The strange 
quark mass is taken as m s (lGeV) = 0.15 GeV, and the cut-off parameter is chosen as 
fj, = 1 GeV. 

Let us evaluate numerically the sum rules for ((y) and (i(y)- The continuum thresholds 
ujq and U\ in (J2~%j) and (}2T?j) are determined by requiring stability of these sum rules. One 



finds that 1.7 GeV < u < 2.2 GeV and 2.6 GeV < u x < 3.1 GeV [11|, |33j, |38j . Imposing 
usual criterion on the ratio of contribution of the higher-order power corrections and that 
of the continuum, we find that for the central values of the condensates given in (|30|). if 
the threshold parameter uo c lies in the range 2.3 < uo c < 2.7 GeV, there is an acceptable 
"stability window" T = 0.8 — 1.2 GeV in which the calculation results do not change 
appreciably. 

The values of the form factors ((y) and Ci(y) at zero recoil as functions of the Borel 
parameter are shown in Fig. 1(a) and 2(a), for three different values of the continuum 
threshold uj c . One can see that the variation is quit moderate in the range 0.8 < T < 1.2 
GeV. The numerical results for ((y) and Ci{y) are shown in Fig. 1(b) and 2(b), where the 
curves refer to three different values of u c and T is fixed at T = 1.0 GeV. 
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0.8 0.9 1 

TGeV 



1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 



FIG. 1: Fig. 1(a) shows the dependence of C(l) on the Borel parameter T for the continuum 
threshold in the range 2.3 < uo c < 2.7 GeV; Fig. 1(b) shows Isgur-Wise function £(y) with T = 1.0 
GeV. 

In order to evaluate numerically the sum rules for ri^ e (y) and t]^ e (y), we need to specif 
the following HQET parameters as input, which are obtained by QCD sum rules 



11 



eciiy 

m 



K. 



3/2 = -(1.6 ± 0.30) GeV 2 
,i /2 = -(1.2 ±0.20) GeV 2 



G 
G 



K 



',3/2 
,1/2 



K 



1.0 ±0.45) GeV 
1.6 ±0.6) GeV. 



(31) 



In Fig. 3(a), the sum rule for rj^ e (y) is plotted at zero recoil as a function of Borel parameter 
for various choices of the continuum thresholds in the range 2.5 < uj c < 2.9. Fig. 3(b) 
shows the y dependence of the form factor rj^ e (y) for the central value of HQET parameters 
and T = 1.2 GeV. It should be noted that apart from the uncertainty from the sum rule 
working window, there is uncertainty to a large extent due to the variation of K's and G' K s 
in the numerical analysis. The numerical evaluation for the sum rule of i]^ e (y) can follow 
the same procedure. 

The numerical analysis shows that all Isgur-Wise functions ((y), (i(y) and rj^ e (y) are 
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FIG. 2: Numerical result for the sum rule (|23|) : (a) dependence of Ci(l) on the Borel parameter 
T for the continuum threshold in the range 2.3 < u c < 2.7 GeV; (b) The form factor Ci(y) with 
T = 1.0 GeV. 
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FIG. 3: Results of the numerical evaluation for the sum rules l|27|l. The curves refer to choices of 
the threshold parameters: uj c = 2.5, tu c = 2.7, tu c = 2.9, from top to bottom. 



slowly varying functions in the allowed kinematic range for B s — > D* £u and B s — > D* sX tv 
decays. They can be well fitted by the linear approximation 

C(y) = C(l) (1 - 0.5(y - 1)) , C(l) = 0.45 ± 0.05 GeV 

Cx(y) = -Ci(l)(l + 0.4(y-l)) , C(l) = 0.65 ±0.06 GeV 

Vl(y) = -Vl(l) (l-0.9(y-l)), i^(l) = 1.7 ±0.2 GeV, 

vL(y) = -?7L(l)(l-0.9(y-l)), 7^,(1) = 1.6 ±0.2 GeV. (32) 

The errors reflect the uncertainty due to co's and T. The uncertainty due to the variation 
of the QCD and HQET parameters is not included, which may reach 5% or more. The 
systematic error resulted from the use of quark-hadron duality above u> c is difficult to 
estimate. Conservatively speaking, there is a 10% systematic error. 

Above parameterizations of the Isgur-Wise functions can be used to calculate the total 
semileptonic rates and decay branching ratios by integrating Eqs. (J!Zaj) and (|2bj) . The 
values of A and A' can be obtained from two-point sum rules (}28|) and (|29|) . respectively. 
They are: A = 0.62 GeV j^] and A' = 0.86 GeV jll|. The quark masses are taken to be 
mj = 4.7 GeV, m c = 1.4 GeV. We use the physical masses, mg = 5.369 |41l, mo. = 2.317 
and Tn £)» = 2.457 [3J], for B s , D* and D* sl mesons. The maximal values of y in the present 
case are y 0max = (1 + rg)/2r = 1.374 and y lmax = (1 + rf)/2n = 1.321. 

In Table U we present our results for decay rates and branching ratios, as well as those 
in the infinitely heavy quark limit. We have taken tb s = 1-46 ps 41]. In the calculation, 
the central values for the Isgur-Wise functions in ()32j) are taken, and the theoretical un- 
certainties are not included. The ratios of the two semileptonic rates for B s decays into 
D% and D *i mesons both in taking account of the I/tuq corrections and in the infinitely 
heavy quark mass limit are 

= B(B S -> D*J,v) _ f 2.05 with l/m Q , 

~ B(B S — > D° sl Iv) \ 1.16 m Q ^oc. 1 ' 
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The branching ratio for B s — > D* £v decay exceeds the one for B s — > D* sl lv in both cases. 







With 1/rriQ 




— » oo 






Decay 


T Br 


r 


Br 


R 


B 




0.38 0.10 


0.31 


0.08 


1.23 


B 


- D* s0 £u 


0.78 0.20 


0.36 


0.09 


2.17 



TABLE I: Decay rates T (in units of \V cb /0. 04| 2 x 10~ 15 GeV) and branching ratios BR (in %) 
for B s — > D^s*^iv decays in taking account of the 1/rriQ corrections and in the infinitely heavy 
quark mass limit. R is a ratio of branching ratios including C(l/mg) corrections to branching 
ratios in the infinitely heavy quark mass limit. 

The numerical predictions in Table |l] indicate that a substantial part of the inclusive 
semileptonic B s decays should go to excited D s meson states. In future hadron B factories 
the D s resonant states can be produced directly in a considerable amount of branching 
ratio from the weak decay of the B s meson. The study of the semileptonic B s decays to 
excited D s states can provide some information about the structure and the properties for 
the newly observed D* sJ states. 

From Table IJ we see that the B — > D* £u decay rate receives large 1/rriQ contributions 
and gets a sharp increase, while the B — > B* sX iv decay rate is only moderately increased 
by subleading 1/rriQ corrections. The reason for this is as following. From Eqs. |Q) and 
(II 1|) we see that the decay matrix elements at zero recoil are determined by form factors 
g+(l) and (7vj.(l)> which receive non- vanishing contributions from first order heavy quark 
mass corrections. Explicitly, 

9+0) = -^c + £6)(A'-A)C(l), 

g Vl {\) = (e c -3e 6 )(A'-A)C(l) . (34) 

At zero recoil the form factor gy 1 is suppressed by a very small factor e c — 3eb ~ 0.04GeV _1 . 
As a result the B s — ► D* sl ev decay rate is only slightly increased by subleading I/tuq 
corrections. On the other hand, B s — > D* eu decay rate receives a large enhancement 
from 1/itiq corrections. Note that the sharp increase of B s — > D* sQ ev decay rate by 1/ttiq 
corrections does not imply the breakdown of the heavy quark expansion. This is because 
the allowed kinematic ranges for B — > Ds$iv is fairly small, the contribution to the decay 
rate of the rather small I/tuq corrections is substantially increased. Hence it is rather a 
result of kinematical and dynamical effects. 

In summary, we have presented the investigation for semileptonic B s decays into excited 
D s mesons. Within the framework of HQET we have applied the QCD sum rules to 
calculate the universal Isgur-Wise functions up to the subleading order of the heavy quark 
expansion. The differential decay widths and the branching ratios for the decays B s — > 
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D* £u and B s — > D* sl lv are computed up to the order of 1/mg corrections. The decay 
rates are substantially influenced by the inclusion of the first order l/rug corrections. 

With the assumption that the newly discovered states _D*j(2317) and _D*j(2460) can be 
identified with spin symmetry doublet (D^D*^, it is worthwhile to study the semileptonic 
B s decays to these newly observed states. In hadron B factories these p-wave excited 
D s states can be produced directly from the semileptonic decay of the B s meson with a 
considerable amount of branching ratio. A measurement of the B s — ► D* sJ can provide 
some information on the nature of the new D* , mesons. 
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